Flexible TAM requirement of TnpB enables efficient single-nucleotide editing with expanded targeting scope.

TnpBs encoded by the IS200/IS605 family transposon are among the most abundant prokaryotic proteins from which type V CRISPR-Cas nucleases may have evolved. Since bacterial TnpBs can be programmed for RNA-guided dsDNA cleavage in the presence of a transposon-adjacent motif (TAM), these nucleases hold immense promise for genome editing. However, the activity and targeting specificity of TnpB in homology-directed gene editing remain unknown. Here we report that a thermophilic archaeal TnpB enables efficient gene editing in the natural host. Interestingly, the TnpB has different TAM requirements for eliciting cell death and for facilitating gene editing. By systematically characterizing TAM variants, we reveal that the TnpB recognizes a broad range of TAM sequences for gene editing including those that do not elicit apparent cell death. Importantly, TnpB shows a very high targeting specificity on targets flanked by a weak TAM. Taking advantage of this feature, we successfully leverage TnpB for efficient single-nucleotide editing with templated repair. The use of different weak TAM sequences not only facilitates more flexible gene editing with increased cell survival, but also greatly expands targeting scopes, and this strategy is probably applicable to diverse CRISPR-Cas systems.

Knockdown of PGBD5 inhibits the malignant progression of glioma through upregulation of the PPAR pathway.

Glioma is the most common type of primary intracranial malignant tumor, and because of its high invasiveness and recurrence, its prognosis remains poor. The present study investigated the biological function of piggyBac transportable element derived 5 (PGBD5) in glioma. Glioma and para-cancerous tissues were obtained from five patients. Reverse transcription-quantitative PCR and western blotting were used to detect the expression levels of PGBD5. Transwell assay and flow cytometry were used to evaluate cell migration, invasion, apoptosis and cell cycle distribution. In addition, a nude mouse tumor transplantation model was established to study the downstream pathways of PGBD5 and the molecular mechanism was analyzed using transcriptome sequencing. The mRNA and protein expression levels of PGBD5 were increased in glioma tissues and cells. Notably, knockdown of PGBD5 in vitro could inhibit the migration and invasion of glioma cells. In addition, the knockdown of PGBD5 expression promoted apoptosis and caused cell cycle arrest in the G2/M phase, thus inhibiting cell proliferation. Furthermore, in vivo experiments revealed that knockdown of PGBD5 expression could inhibit Ki67 expression and slow tumor growth. Changes in PGBD5 expression were also shown to be closely related to the peroxisome proliferator-activated receptor (PPAR) signaling pathway. In conclusion, interference with PGBD5 could inhibit the malignant progression of glioma through the PPAR pathway, suggesting that PGBD5 may be a potential molecular target of glioma.

CRISPR RNA-Guided Transposases Facilitate Dispensable Gene Study in Phage.

Phages provide a potential therapy for multi-drug-resistant (MDR) bacteria. However, a significant portion of viral genes often remains unknown, posing potential dangers. The identification of non-essential genes helps dissect and simplify phage genomes, but current methods have various limitations. In this study, we present an in vivo two-plasmid transposon insertion system to assess the importance of phage genes, which is based on the V. cholerae transposon Tn6677, encoding a nuclease-deficient type I-F CRISPR-Cas system. We first validated the system in Pseudomonas aeruginosa PAO1 and its phage S1. We then used the selection marker AcrVA1 to protect transposon-inserted phages from CRISPR-Cas12a and enriched the transposon-inserted phages. For a pool of selected 10 open-reading frames (2 known functional protein genes and 8 hypothetical protein genes) of phage S1, we identified 5 (2 known functional protein genes and 3 hypothetical protein genes) as indispensable genes and the remaining 5 (all hypothetical protein genes) as dispensable genes. This approach offers a convenient, site-specific method that does not depend on homologous arms and double-strand breaks (DSBs), holding promise for future applications across a broader range of phages and facilitating the identification of the importance of phage genes and the insertion of genetic cargos.

Iterative assay for transposase-accessible chromatin by sequencing to isolate functionally relevant neuronal subtypes.

The Drosophila brain contains tens of thousands of distinct cell types. Thousands of different transgenic lines reproducibly target specific neuron subsets, yet most still express in several cell types. Furthermore, most lines were developed without a priori knowledge of where the transgenes would be expressed. To aid in the development of cell type-specific tools for neuronal identification and manipulation, we developed an iterative assay for transposase-accessible chromatin (ATAC) approach. Open chromatin regions (OCRs) enriched in neurons, compared to whole bodies, drove transgene expression preferentially in subsets of neurons. A second round of ATAC-seq from these specific neuron subsets revealed additional enriched OCR2s that further restricted transgene expression within the chosen neuron subset. This approach allows for continued refinement of transgene expression, and we used it to identify neurons relevant for sleep behavior. Furthermore, this approach is widely applicable to other cell types and to other organisms.

Childhood cancer mutagenesis caused by transposase-derived PGBD5.

Genomic rearrangements are a hallmark of most childhood tumors, including medulloblastoma, one of the most common brain tumors in children, but their causes remain largely unknown. Here, we show that PiggyBac transposable element derived 5 (Pgbd5) promotes tumor development in multiple developmentally accurate mouse models of Sonic Hedgehog (SHH) medulloblastoma. Most Pgbd5-deficient mice do not develop tumors, while maintaining normal cerebellar development. Ectopic activation of SHH signaling is sufficient to enforce cerebellar granule cell progenitor-like cell states, which exhibit Pgbd5-dependent expression of distinct DNA repair and neurodevelopmental factors. Mouse medulloblastomas expressing Pgbd5 have increased numbers of somatic structural DNA rearrangements, some of which carry PGBD5-specific sequences at their breakpoints. Similar sequence breakpoints recurrently affect somatic DNA rearrangements of known tumor suppressors and oncogenes in medulloblastomas in 329 children. This identifies PGBD5 as a medulloblastoma mutator and provides a genetic mechanism for the generation of oncogenic DNA rearrangements in childhood cancer.

An optimized polymeric delivery system for piggyBac transposition.

The piggyBac transposon/transposase system has been explored for long-term, stable gene expression to execute genomic integration of therapeutic genes, thus emerging as a strong alternative to viral transduction. Most studies with piggyBac transposition have employed physical methods for successful delivery of the necessary components of the piggyBac system into the cells. Very few studies have explored polymeric gene delivery systems. In this short communication, we report an effective delivery system based on low molecular polyethylenimine polymer with lipid substitution (PEI-L) capable of delivering three components, (i) a piggyBac transposon plasmid DNA carrying a gene encoding green fluorescence protein (PB-GFP), (ii) a piggyBac transposase plasmid DNA or mRNA, and (iii) a 2 kDa polyacrylic acid as additive for transfection enhancement, all in a single complex. We demonstrate an optimized formulation for stable GFP expression in two model cell lines, MDA-MB-231 and SUM149 recorded till day 108 (3.5 months) and day 43 (1.4 months), respectively, following a single treatment with very low cell number as starting material. Moreover, the stability of the transgene (GFP) expression mediated by piggyBac/PEI-L transposition was retained following three consecutive cryopreservation cycles. The success of this study highlights the feasibility and potential of employing a polymeric delivery system to obtain piggyBac-based stable expression of therapeutic genes.

Application of Single-Cell Assay for Transposase-Accessible Chromatin with High Throughput Sequencing in Plant Science: Advances, Technical Challenges, and Prospects.

The Single-cell Assay for Transposase-Accessible Chromatin with high throughput sequencing (scATAC-seq) has gained increasing popularity in recent years, allowing for chromatin accessibility to be deciphered and gene regulatory networks (GRNs) to be inferred at single-cell resolution. This cutting-edge technology now enables the genome-wide profiling of chromatin accessibility at the cellular level and the capturing of cell-type-specific cis-regulatory elements (CREs) that are masked by cellular heterogeneity in bulk assays. Additionally, it can also facilitate the identification of rare and new cell types based on differences in chromatin accessibility and the charting of cellular developmental trajectories within lineage-related cell clusters. Due to technical challenges and limitations, the data generated from scATAC-seq exhibit unique features, often characterized by high sparsity and noise, even within the same cell type. To address these challenges, various bioinformatic tools have been developed. Furthermore, the application of scATAC-seq in plant science is still in its infancy, with most research focusing on root tissues and model plant species. In this review, we provide an overview of recent progress in scATAC-seq and its application across various fields. We first conduct scATAC-seq in plant science. Next, we highlight the current challenges of scATAC-seq in plant science and major strategies for cell type annotation. Finally, we outline several future directions to exploit scATAC-seq technologies to address critical challenges in plant science, ranging from plant ENCODE(The Encyclopedia of DNA Elements) project construction to GRN inference, to deepen our understanding of the roles of CREs in plant biology.

Mage transposon: a novel gene delivery system for mammalian cells.

Transposons, as non-viral integration vectors, provide a secure and efficient method for stable gene delivery. In this study, we have discovered Mage (MG), a novel member of the piggyBac(PB) family, which exhibits strong transposability in a variety of mammalian cells and primary T cells. The wild-type MG showed a weaker insertion preference for near genes, transcription start sites (TSS), CpG islands, and DNaseI hypersensitive sites in comparison to PB, approaching the random insertion pattern. Utilizing in silico virtual screening and feasible combinatorial mutagenesis in vitro, we effectively produced the hyperactive MG transposase (hyMagease). This variant boasts a transposition rate 60% greater than its native counterpart without significantly altering its insertion pattern. Furthermore, we applied the hyMagease to efficiently deliver chimeric antigen receptor (CAR) into T cells, leading to stable high-level expression and inducing significant anti-tumor effects both in vitro and in xenograft mice models. These findings provide a compelling tool for gene transfer research, emphasizing its potential and prospects in the domains of genetic engineering and gene therapy.

Increasing transposase abundance with ocean depth correlates with a particle-associated lifestyle.

Transposases are mobile genetic elements that move within and between genomes, promoting genomic plasticity in microorganisms. In marine microbial communities, the abundance of transposases increases with depth, but the reasons behind this trend remain unclear. Our analysis of metagenomes from the Tara Oceans and Malaspina Expeditions suggests that a particle-associated lifestyle is the main covariate for the high occurrence of transposases in the deep ocean, and this trend holds true for individual genomes as well as in a community-wide sense. We observed a strong and depth-independent correlation between transposase abundance and the presence of biofilm-associated genes, as well as the prevalence of secretory enzymes. This suggests that mobile genetic elements readily propagate among microbial communities within crowded biofilms. Furthermore, we show that particle association positively correlates with larger genome size, which is in turn associated with higher transposase abundance. Cassette sequences associated with transposons are enriched with genes related to defense mechanisms, which are more highly expressed in the deep sea. Thus, while transposons spread at the expense of their microbial hosts, they also introduce novel genes and potentially benefit the hosts in helping to compete for limited resources. Overall, our results suggest a new understanding of deep ocean particles as highways for gene sharing among defensively oriented microbial genomes.IMPORTANCEGenes can move within and between microbial genomes via mobile genetic elements, which include transposases and transposons. In the oceans, there is a puzzling increase in transposase abundance in microbial genomes as depth increases. To gain insight into this trend, we conducted an extensive analysis of marine microbial metagenomes and metatranscriptomes. We found a significant correlation between transposase abundance and a particle-associated lifestyle among marine microbes at both the metagenome and genome-resolved levels. We also observed a link between transposase abundance and genes related to defense mechanisms. These results suggest that as microbes become densely packed into crowded particles, mobile genes are more likely to spread and carry genetic material that provides a competitive advantage in crowded habitats. This may enable deep sea microbes to effectively compete in such environments.

Application of RNA interference and protein localization to investigate housekeeping and developmentally regulated genes in the emerging model protozoan Paramecium caudatum.

Unicellular eukaryotes represent tremendous evolutionary diversity. However, the molecular mechanisms underlying this diversity remain largely unexplored, partly due to a limitation of genetic tools to only a few model species. Paramecium caudatum is a well-known unicellular eukaryote with an unexpectedly large germline genome, of which only two percent is retained in the somatic genome following sexual processes, revealing extensive DNA elimination. However, further progress in understanding the molecular mechanisms governing this process is hampered by a lack of suitable genetic tools. Here, we report the successful application of gene knockdown and protein localization methods to interrogate the function of both housekeeping and developmentally regulated genes in P. caudatum. Using these methods, we achieved the expected phenotypes upon RNAi by feeding, and determined the localization of these proteins by microinjection of fusion constructs containing fluorescent protein or antibody tags. Lastly, we used these methods to reveal that P. caudatum PiggyMac, a domesticated piggyBac transposase, is essential for sexual development, and is likely to be an active transposase directly involved in DNA cleavage. The application of these methods lays the groundwork for future studies of gene function in P. caudatum and can be used to answer important biological questions in the future.

Transposition mechanism of ISApl1-the determinant of colistin resistance dissemination.

Multidrug-resistant Enterobacteriaceae, a prominent family of gram-negative pathogenic bacteria, causes a wide range of severe diseases. Strains carrying the mobile colistin resistance (mcr-1) gene show resistance to polymyxin, the last line of defense against multidrug-resistant gram-negative bacteria. However, the transmission of mcr-1 is not well understood. In this study, genomes of mcr-1-positive strains were obtained from the NCBI database, revealing their widespread distribution in China. We also showed that ISApl1, a crucial factor in mcr-1 transmission, is capable of self-transposition. Moreover, the self-cyclization of ISApl1 is mediated by its own encoded transposase. The electrophoretic mobility shift assay experiment validated that the transposase can bind to the inverted repeats (IRs) on both ends, facilitating the cyclization of ISApl1. Through knockout or shortening of IRs at both ends of ISApl1, we demonstrated that the cyclization of ISApl1 is dependent on the sequences of the IRs at both ends. Simultaneously, altering the ATCG content of the bases at both ends of ISApl1 can impact the excision rate by modifying the binding ability between IRs and ISAPL1. Finally, we showed that heat-unstable nucleoid protein (HU) can inhibit ISApl1 transposition by binding to the IRs and preventing ISAPL1 binding and expression. In conclusion, the regulation of ISApl1-self-circling is predominantly controlled by the inverted repeat (IR) sequence and the HU protein. This molecular mechanism deepens our comprehension of mcr-1 dissemination.

Bacterial genome engineering using CRISPR-associated transposases.

Clustered regularly interspaced short palindromic repeats (CRISPR)-associated transposases have the potential to transform the technology landscape for kilobase-scale genome engineering, by virtue of their ability to integrate large genetic payloads with high accuracy, easy programmability and no requirement for homologous recombination machinery. These transposons encode efficient, CRISPR RNA-guided transposases that execute genomic insertions in Escherichia coli at efficiencies approaching ~100%. Moreover, they generate multiplexed edits when programmed with multiple guides, and function robustly in diverse Gram-negative bacterial species. Here we present a detailed protocol for engineering bacterial genomes using CRISPR-associated transposase (CAST) systems, including guidelines on the available vectors, customization of guide RNAs and DNA payloads, selection of common delivery methods, and genotypic analysis of integration events. We further describe a computational CRISPR RNA design algorithm to avoid potential off-targets, and a CRISPR array cloning pipeline for performing multiplexed DNA insertions. The method presented here allows the isolation of clonal strains containing a novel genomic integration event of interest within 1-2 weeks using available plasmid constructs and standard molecular biology techniques.

Targeted DNA integration in human cells without double-strand breaks using CRISPR-associated transposases.

Conventional genome engineering with CRISPR-Cas9 creates double-strand breaks (DSBs) that lead to undesirable byproducts and reduce product purity. Here we report an approach for programmable integration of large DNA sequences in human cells that avoids the generation of DSBs by using Type I-F CRISPR-associated transposases (CASTs). We optimized DNA targeting by the QCascade complex through protein design and developed potent transcriptional activators by exploiting the multi-valent recruitment of the AAA+ ATPase TnsC to genomic sites targeted by QCascade. After initial detection of plasmid-based integration, we screened 15 additional CAST systems from a wide range of bacterial hosts, identified a homolog from Pseudoalteromonas that exhibits improved activity and further increased integration efficiencies. Finally, we discovered that bacterial ClpX enhances genomic integration by multiple orders of magnitude, likely by promoting active disassembly of the post-integration CAST complex, akin to its known role in Mu transposition. Our work highlights the ability to reconstitute complex, multi-component machineries in human cells and establishes a strong foundation to exploit CRISPR-associated transposases for eukaryotic genome engineering.

AAV-mediated delivery of a Sleeping Beauty transposon and an mRNA-encoded transposase for the engineering of therapeutic immune cells.

Engineering cells for adoptive therapy requires overcoming limitations in cell viability and, in the efficiency of transgene delivery, the duration of transgene expression and the stability of genomic integration. Here we report a gene-delivery system consisting of a Sleeping Beauty (SB) transposase encoded into a messenger RNA delivered by an adeno-associated virus (AAV) encoding an SB transposon that includes the desired transgene, for mediating the permanent integration of the transgene. Compared with lentiviral vectors and with the electroporation of plasmids of transposon DNA or minicircle DNA, the gene-delivery system, which we named MAJESTIC (for 'mRNA AAV-SB joint engineering of stable therapeutic immune cells'), offers prolonged transgene expression, as well as higher transgene expression, therapeutic-cell yield and cell viability. MAJESTIC can deliver chimeric antigen receptors (CARs) into T cells (which we show lead to strong anti-tumour activity in vivo) and also transduce natural killer cells, myeloid cells and induced pluripotent stem cells with bi-specific CARs, kill-switch CARs and synthetic T-cell receptors.

Insertion sequence excision is enhanced by a protein that catalyzes branch migration and promotes microhomology-mediated end joining.

Insertion sequence (IS)-excision enhancer (IEE) promotes the excision of ISs in the genome of enterohemorrhagic Escherichia coli O157. Because IEE-dependent IS excision occurs in the presence of transposase, the process of IS transposition may be involved in IS excision; however, little is understood about the molecular mechanisms of IS excision. Our in vitro analysis revealed that IEE exhibits DNA-dependent ATPase activity, which is activated by branched DNA. IEE also catalyzes the branch migration of fork-structured DNA. These results suggest that IEE remodels branched structures of the IS transposition intermediate. Sequence analysis of recombination sites in IS-excision products suggested that microhomologous sequences near the ends of the IS are involved in IS excision. IEE promoted microhomology-mediated end joining (MMEJ), in which base pairing between 6-nucleotides complementary ends of two 3'-protruding DNAs and subsequent elongation of the paired DNA strand occurred. IS-excision frequencies were significantly decreased in cells producing IEE mutants that had lost either branch migration or MMEJ activity, which suggests that these activities of IEE are required for IS excision. Based on our results, we propose a model for IS excision triggered by IEE and transposase.

Chromosome recombination and modification by LoxP-mediated evolution in Vibrio natriegens using CRISPR-associated transposases.

Chromosome rearrangement by LoxP-mediated evolution has emerged as a powerful approach to studying how chromosome architecture impacts phenotypes. However, it relies on the in vitro synthesis of artificial chromosomes. The recently reported CRISPR-associated transposases (CASTs) held great promise for the efficient insertion of abundant LoxP sites directly onto the genome of wild-type strains. In this study, with the fastest-growing bacterium Vibrio natrigens (V. natriegens) as an object, a multiplex genome integration tool derived from CASTs was employed to achieve the insertion of cargo genes at eight specific genomic loci within 2 days. Next, we introduced 30 LoxP sites onto chromosome 2 (Chr2) of V. natriegens. Rigorously induced Cre recombinase was used to demonstrate Chromosome Rearrangement and Modification by LoxP-mediated Evolution (CRaMbLE). Growth characterization and genome sequencing showed that the ~358 kb fragment on Chr2 was accountable for the rapid growth of V. natriegens. The enabling tools we developed can help identify genomic regions that influence the rapid growth of V. natriegens without a prior understanding of genome mechanisms. This groundbreaking demonstration may also be extended to other organisms such as Escherichia coli, Pseudomonas putida, Bacillus subtilis, and so on.

The horizontal gene transfer of perchlorate reduction genomic island in three bacteria from an ecological niche.

Three new strains of dissimilatory perchlorate-reducing bacteria (DPRB), QD19-16, QD1-5, and P3-1, were isolated from an active sludge. Phylogenetic trees based on 16S rRNA genes indicated that QD19-16, QD1-5, and P3-1 belonged to Brucella, Acidovorax, and Citrobacter, respectively, expanding the distribution of DPRB in the Proteobacteria. The three strains were gram-negative and facultative anaerobes with rod-shaped cells without flagella, which were 1.0-1.6 µm long and 0.5-0.6 µm wide. The three DPRB strains utilized similar broad spectrum of electron donors and acceptors and demonstrated a similar capability to reduce perchlorate within 6 days. The enzyme activity of perchlorate reductase in QD19-16 toward chlorate was higher than that toward perchlorate. The high sequence similarity of the perchlorate reductase operon and chlorite dismutase genes in the perchlorate reduction genomic islands (PRI) of the three strains implied that they were monophyletic origin from a common ancestral PRI. Two transposase genes (tnp1 and tnp2) were found in the PRIs of strain QD19-16 and QD1-5, but were absent in the strain P3-1 PRI. The presence of fragments of IR sequences in the P3-1 PRI suggested that P3-1 PRI had previously contained these two tnp genes. Therefore, it is plausible to suggest that a common ancestral PRI transferred across the strains Brucella sp. QD19-16, Acidovorax sp. QD1-5, and Citrobacter sp. P3-1 through horizontal gene transfer, facilitated by transposases. These results provided a direct evidence of horizontal gene transfer of PRI that could jump across phylogenetically unrelated bacteria through transposase. KEY POINTS: • Three new DPRB strains can effectively remove high concentration of perchlorate. • The PRIs of three DPRB strains are acquired from a single ancestral PRI. • PRIs are incorporated into different bacteria genome through HGT by transposase.

Dynamics of chromatin accessibility and genome wide control of desiccation tolerance in the resurrection plant Haberlea rhodopensis.

BACKGROUND: Drought is one of the main consequences of global climate change and this problem is expected to intensify in the future. Resurrection plants evolved the ability to withstand the negative impact of long periods of almost complete desiccation and to recover at rewatering. In this respect, many physiological, transcriptomic, proteomic and genomic investigations have been performed in recent years, however, few epigenetic control studies have been performed on these valuable desiccation-tolerant plants so far. RESULTS: In the present study, for the first time for resurrection plants we provide evidences about the differential chromatin accessibility of Haberlea rhodopensis during desiccation stress by ATAC-seq (Assay for Transposase Accessible Chromatin with high-throughput sequencing). Based on gene similarity between species, we used the available genome of the closely related resurrection plant Dorcoceras hygrometricum to identify approximately nine hundred transposase hypersensitive sites (THSs) in H. rhodopensis. The majority of them corresponds to proximal and distal regulatory elements of different genes involved in photosynthesis, carbon metabolism, synthesis of secondary metabolites, cell signalling and transcriptional regulation, cell growth, cell wall, stomata conditioning, chaperons, oxidative stress, autophagy and others. Various types of binding motifs recognized by several families of transcription factors have been enriched from the THSs found in different stages of drought. Further, we used the previously published RNA-seq data from H. rhodopensis to evaluate the expression of transcription factors putatively interacting with the enriched motifs, and the potential correlation between the identified THS and the expression of their corresponding genes. CONCLUSIONS: These results provide a blueprint for investigating the epigenetic regulation of desiccation tolerance in resurrection plant H. rhodopensis and comparative genomics between resurrection and non-resurrection species with available genome information.

Mechanism of target site selection by type V-K CRISPR-associated transposases.

CRISPR-associated transposases (CASTs) repurpose nuclease-deficient CRISPR effectors to catalyze RNA-guided transposition of large genetic payloads. Type V-K CASTs offer potential technology advantages but lack accuracy, and the molecular basis for this drawback has remained elusive. Here, we reveal that type V-K CASTs maintain an RNA-independent, "untargeted" transposition pathway alongside RNA-dependent integration, driven by the local availability of TnsC filaments. Using cryo-electron microscopy, single-molecule experiments, and high-throughput sequencing, we found that a minimal, CRISPR-less transpososome preferentially directs untargeted integration at AT-rich sites, with additional local specificity imparted by TnsB. By exploiting this knowledge, we suppressed untargeted transposition and increased type V-K CAST specificity up to 98.1% in cells without compromising on-target integration efficiency. These findings will inform further engineering of CAST systems for accurate, kilobase-scale genome engineering applications.

Improved quality metrics for association and reproducibility in chromatin accessibility data using mutual information.

BACKGROUND: Correlation metrics are widely utilized in genomics analysis and often implemented with little regard to assumptions of normality, homoscedasticity, and independence of values. This is especially true when comparing values between replicated sequencing experiments that probe chromatin accessibility, such as assays for transposase-accessible chromatin via sequencing (ATAC-seq). Such data can possess several regions across the human genome with little to no sequencing depth and are thus non-normal with a large portion of zero values. Despite distributed use in the epigenomics field, few studies have evaluated and benchmarked how correlation and association statistics behave across ATAC-seq experiments with known differences or the effects of removing specific outliers from the data. Here, we developed a computational simulation of ATAC-seq data to elucidate the behavior of correlation statistics and to compare their accuracy under set conditions of reproducibility. RESULTS: Using these simulations, we monitored the behavior of several correlation statistics, including the Pearson's R and Spearman's [Formula: see text] coefficients as well as Kendall's [Formula: see text] and Top-Down correlation. We also test the behavior of association measures, including the coefficient of determination R[Formula: see text], Kendall's W, and normalized mutual information. Our experiments reveal an insensitivity of most statistics, including Spearman's [Formula: see text], Kendall's [Formula: see text], and Kendall's W, to increasing differences between simulated ATAC-seq replicates. The removal of co-zeros (regions lacking mapped sequenced reads) between simulated experiments greatly improves the estimates of correlation and association. After removing co-zeros, the R[Formula: see text] coefficient and normalized mutual information display the best performance, having a closer one-to-one relationship with the known portion of shared, enhanced loci between simulated replicates. When comparing values between experimental ATAC-seq data using a random forest model, mutual information best predicts ATAC-seq replicate relationships. CONCLUSIONS: Collectively, this study demonstrates how measures of correlation and association can behave in epigenomics experiments. We provide improved strategies for quantifying relationships in these increasingly prevalent and important chromatin accessibility assays.